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Divide Geo-operator Volume 
Into N Cross-Sections 



Threshold Regions as Three Materials 

1 ) exterior physical matrix 

2) natural resource 

3) nonessential material 



Apply Onion Skin Algorithm 

find layers of each region to find contiguity by rings 
find max and min rings of each material 



I 



Determine: 

Do the max physical matrix ring 
and the min natural resource ring intersect? 
and 

Does the max natural resource ring and any 
nonessential material ring make contact? 



No 

^assign val ue of 1 



assign value of 0 




^Sum 

values for all 
cross-sections 
and divide final 
result by N 




Figure 18 
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ProgramName Follower 18/21 
Begin 

^Commenting Select the path 
2100 Input N, pathpoints 

2110 Set startpoint < Pathpoints(l) 

2120 Set endpoint < Pathpoints(N) 
2130 Set Refpoint < startpoint 

2140 Set voxeloistance < pathpoints(2)- Pathpoints(l) . 

^Commenting Select the geo-operator 
2200 Input xiength, ylength, zlength 
2210 Input voxref^ 
2220 input Direction 
2230 input Algorithm 
2240 Input validitylnfringementRule 
2250 Input K, DecisionBoundaryTable 
2260 Set opx< xiength 
2270 Set opy< ylength 
2280 Set opz< zlength 

2300 LoopSet M < N ; Dowhile M > 0 

^Commenting Move the aeo-ooeratorif not at the startpoint 
2310 If (Refpoint Equal startpoint) 

2320 Set Newpoint < Refpoint 

Otherwise 

2330 Set Newpoint < Refpoint + voxelDistance 

EndIf 

2340 Calculate smoothRotatedmdices(pathpoints opx,opy,opz, ,Direction,RefpointNewpoint) 
2350 Calculate verticalFloat(pathpoints J opx J opy l opz 1 Direction > RefpointNewpoint) 

^Commenting Determine if valid computation is possible 
2400 Calculate validityinfringementRule(SmoothRotatedlndices, VerticalFloat, Op xyz Newpoint) 
2410 If (validitylnfringementRule <= 0) 

2420 Set voxelvalidityAttribute(pathpoints(M)) < False 

Otherwise 

2430 Set voxeivalidityAttribute(pathpoints(M)) < True 

2440 LoopSet H < size(opx,opy,opz) ; Dowhile H > 0 

^Commenting Calculate and store the aeo-operator 
2500 Input Brickofbata(smoothRotatedmdices) 

2510 Calculate Algorithm(Pathpoint(M)) 

2520 Store calculationResult(M) 

Calculate DedsionBoundai7Result(DecisionBoundaryTable(calculationResult(M))) 
f^Z Set voxelPathAttribute(Pathpoints(M)) < Pathcolor 

2540 Set voxelcalculationAttribute(Pathpoints(M)) ( 

2560 DecisionBoundaryResult(M) 
2570 Store voxel( . . . )Attribute(pathpoints(M)) 

2580 Decrement H 

2590 EndofLoopOn H 

EndIf 

^Commenting Increment the position and ao to next path point 
2600 Set Refpoint < Refpoint + voxelDistance 
2610 Decrement M 
2620 EndofLoopOn M 

-Commenting Visualize results 
2700 Recall voxel(...)Attribute(pathpoint(...)) 

2710 Display voxel(...)Attribute(pathpoint(...)j r-r^...^^ *%A 

endprogram. ilQUrQ jci 
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^ current position^ 



Record feature 
position 
i zz 
Record feature 

absolute size 



Choose next voxel 
candidate along path 

Select all data 
needed to cover 
operator extents 

i Z 




Perform 

feature 

correlation 



Record vertical 
translation of centroid 



Record vertical and 
horizontal scaling 



Loop until 




operator 




^ complete ^ 
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Loop until 




path 




^complete^ 
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Figure 24 



Apply Geo-operator 
to Path 




